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ABSTRACT 

Observations of magnetism in very low mass stars recently made important progress, revealing 
characteristics that are now to be understood in the framework of dynamo theory. In parallel, 
there is growing evidence that dynamo processes in these stars share many similarities with 
planetary dynamos. We investigate the extent to which the weak vs strong field bistability 
predicted for the geodynamo can apply to recent observations of two groups of very low 
mass fully-convective stars sharing similar stellar parameters but generating radically different 
types of magnetic fields. Our analysis is based on previously published spectropolarimetric 
and spectroscopic data. We argue that these can be interpreted in the framework of weak and 
strong field dynamos. 

Key words: Dynamo -Planets and satellites: magnetic fields -Stars: magnetic fields -Stars: 
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1 INTRODUCTION 

Many stars possessing an outer convective envelope - like the Sun 
- exhibit a variety of activity phenomena (e.g., cool spots produc- 
ing photometric variations, a hot corona detected at radio and X-ray 
wavelengths) powered by their magnetic field. The latter is gener- 
ated by dynamo processes converting kinetic energy (due to ther- 
mal convection) into magnetic energy. In the Sun, and other solar- 
type stars, the tachocline, a thin shear layer at the base of the so- 
lar convection zone, is thought to play an important role in gener - 
ating the magnetic field (e.g., [charbonneau & MacG regoill 19971) . 
On the contrary, main-sequence stars below 0.35 M Q being fully- 
convective, do not possess a tachocline. Dynamo processes in these 
objects are therefore believed to differ significantly from those in 
the Sun; in partic ular, they may operate throughout the wh ole stel- 
lar interior (e.g., Ichabrier & KukeJl2006l ; [Browninsll2008h . 

Measurements of surface magnetic fields on a number of M 
dwarfs (0.08 < M t /M < 0.7) have recently been available us- 
ing two complementary approaches. One is based on spectroscopy: 
the average value of the scalar magnetic field at the surface of 
the star is inferr ed from the analysis of the Zeem an broadening 
of spectral lines JSaaill 19881 : iReiners & Basrill2006l) . The other ap- 
proach uses time-series of circularly polarised spectra and tomo- 
graphic imaging techniques to produce spatially resolved maps 
of the large-scale component (typically up to spherical harmonic 
degree £ msa in the range 6-30, depending on the rotational ve- 
locity) of the vector magnetic field (Zeeman-Doppler Imaging 
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ISemellll989i; iDonati & Browrj|l997l : IDonati et alj |"2006b). We re- 
fer to lMorin et ail ( l2010h and references therein for a more detailed 
comparison. The first spectropol arimetric observations of a fully- 
convec tive M dwarf (V374 Peg, IDonati et al]|2006al : iMorin et ail 
2008a) have revealed that these objects can host magnetic fields 
featuring a long-lived strong dipolar component almost aligned 
with the rotation axis, much more reminiscent of planetary mag- 
netic fields than th ose observed on the Sun or solar-type stars (e.g., 
IDonati et all2003l) . 

Despite many differences between planetary and stellar interi- 
ors, a few recent studies have strengthened the idea that some fun- 
damental properties of stellar dynamos could be captured by sim- 
ple Boussinesq models (i.e. without taking into account the radial 
dependency of the stellar density). An accurate description of the 
interior dynamics of stars requires considering their density strati- 
fication. This can be more reliably achieved by anelastic models. It 
seems however, that some characteristics are ro bust enough to be 
alread y captured by a Boussinesq description. iGoudard & Dormvl 
(2008) have shown, for one given set of parameters, that numer- 
ical simulations in the Boussinesq approximation can reproduce 
some basic characteristics of either the magnetism of planets or 
partly-convective stars - steady axial dipole versus cyclic dynamo 
waves respectively - by varying a single parameter: the relative 
width of the conv e ction zone (a thin shell leads to dynamo waves). 
Christensen et al. (2009) (C09) showed that a scaling law for the 
magnetic field strength originally derived from a large number of 
Boussinesq geodynamo simulations is also applicable to rapidly- 
rotating, fully-convective stars (either main-sequence M dwarfs or 
young contracting T Tauri stars). 

iFeatherstone et al. (2009) show that dynamo action in a fully- 
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convective sphere (simulating the core of a A-type star) can be 
strongly enhanced by using a strong initial dipolar field, as in con- 
trast to the case where only a small seed field is initially present. 
The dipolar solution bears strong similarities with Boussinesq geo- 
dynamo models, as noted by the authors. Moreover the existence 
of two co-existing solutions for a given p arameter set seems to be 
reminiscent of the bistability described in ISimitev & Bussej J2009t) 

and found in Boussinesq models. 

Re cent spectropolarimetric observations by Mori n et all 
(2010a) (hereafter MlOa) have revealed two radically different 
types of large-scale magnetism for M dwarfs with similar masses 
and rotation periods. One possible explanation for these observa- 
tions could be the existence of two dynamo branches in this param- 
eter regime. In the present letter, we briefly summarize the obser- 
vational results on fully-convective stars, describe the theoretical 
framework of the weak vs strong field dynamo bistability, and dis- 
cuss its applicability to very low mass stars. 



2 LARGE-SCALE VERSUS TOTAL MAGNETIC FIELD 

Following the study of V374 Peg, a first spectropolarimetric multi- 
epoch survey was initiated for a sample of 23 active M dwarfs. The 
survey intended to constrain the effects of the s hift towards a fully- 
convective interna l structure on dynamo action dMorin e t al. 2008b; 
iDonati et alj2008l , MlOa), by mapping the large-scale component 
of the surface magnetic field, and assessing the corresponding mag- 
netic flux Q (with a typical uncertainty of the order of 20 %). 

All the fully-convective stars of the sample lie in the so called 
saturated dynamo regime - corresponding to P mt < (5, 10} d for 
a {0.3 5,0.15} M G star (see iKiraga & Stepierj|2007l : IReiners et al.l 
2009). In this regime, the rapid dependence of the magnetic flux 
with the rotation rate observed for slower rotators (Rossby num- 
bers larger than 0.1) suddenly stops. The observations (in terms of 
X-ray activity or Zeeman broadening) are consistent with a mag- 
netic field almost independent of the rotation rate. We will come 
back to this point in section [4] We indeed verified that for stars 
for which such measurements exist, the total magnetic flux inferred 
from unpolarised spectra is in the range of 1-4 kG (i.e. matching 
the C09 scaling law). Two radically different types of large-scale 
magnetic fields are observed, either a strong and steady axial dipole 
field (hereafter SD) or a weaker multipolar, non-axisymmetric field 
configuration in rapid evolution (hereafter WM), whereas no dis- 
tinction between these two groups of stars can be made on the basis 
of mass and rotation only, see Fig.Q] 

All stars in the strong dipole regime have a typical large-scale 
magnetic flux of 1 kG (values comprised between 0.5 and 1.6 kG) 
whereas for those in the weak, multipolar regime the typical value 
is 0.1 kG (all values are lower than 0.2 kG) and it is much more 
variable for a given object. The latter is only found in the parameter 
range M t < 0.15 M G and P mt < 1.5 d in our sample, though the 
limits of the domain in which this behaviour occurs are not yet well 
defined (see Fig. 1 1 121 . additional observations on an larger sample 
of very low mass stars are needed to specify this point. 

Measurements of the total magnetic flux BI from unpolarised 
spectroscopy are available for a number of stars in our spectropo- 
larimetric sampl e (see Fig. [3}, with typical uncertainties in the 
range 0.5-1 kG dReiners et alj|2009l ; IReiners & Basrill2010l) . The 

1 Here we ter m "magnetic fl ux" the average modulus of the surface mag- 
netic field, see iReinersI iipioh for a discussion on this term. 
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Figure 1. Mass-period diagram of fully-convective stars derived from spec- 
tropolarimetric data an d Zeeman-Doppler Imaging (ZDI) by citeMorin08b, 
Phan-Bao et al. 1 2009), and MlOa. Symbol size represents the reconstructed 
magnetic energy, the color ranges from blue to red for purely toroidal to 
purely poloidal field, and the shape depicts the degree of axisymmetry from 
a sharp star for non-axisymmetric to a regular decagon for axisymmetric. 
For a few stars of the sample MlOa could not perform a definite ZDI recon- 
struction, in these cases only an upper limit of the rotation period is known 
and the magnetic flux is extrapolated, those objects are depicted as empty 
symbols. The theoretical fully-convective limit is depicted as a horizontal 
dashed line. Thin solid lines represent contours of constant Rossby number 
Ro=0.01 (left) and 0.1 (right), as estimated in MlOa. 

three objects featuring Si ~ 4 kG are in the SD regime - large red 
decagons - but both SD and WM large-scale magnetic fields are 
found among the stars having Sj ~ 2 kG (see Fig. [3}. We there- 
fore conclude that there is no systematic correlation between the 
unpolarised magnetic flux Si and the large-scale magnetic topol- 
ogy inferred from spectropolarimetric observations. Hence, the two 
different types of magnetic field configurations are only detected 
when considering the large-scale component (probed by spectropo- 
larimetry, and which represents 15-30 % of the total flux in the SD 
regime, but only a few percent in the WM regime) and not the total 
magnetic flux derived from unpolarised spectroscopy. 



3 WEAK AND STRONG FIELD DYNAMOS 

It has been known since IChandrasekhan ( ll96lh that both, magnetic 
fields and rotation, taken separately tend to inhibit convective mo- 
tions, but that if both effects are combined the impeding influences 
of the Lorentz and of the Coriolis force may be partly relaxed, al- 
lowing convection to set in at lower R ayleigh number and to de - 
velop on larger length scales (see also lEltaveb & Roberts! fl970T) . 
This mutual counteraction of rotation and magnetism is most effec- 
tive if the Lorentz and the Coriolis forces are of th e same order of 
magnitude, this is the magnetostrophic balance (see Chandr asekhad 
1 196lt ISoward1ll979l : Istevensonll 1979h . This led lRobertsI Jl978h to 
conjecture the existence of two different dynamo regimes - a weak 
and a strong field branch - and that these different dynamo so- 
lutions could co-exist over some range of parameters (see also 
lRobertdfl988l for a review). The anticipated bifurcation diagram 
(adapted from iRobertsll 19881) is presented in Fig. [4] For dynamos 
belonging to the weak field branch, the Proudman-Taylor constraint 
can only be broken owing to the presence of the viscous or the iner- 
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Figure 2. Average large-scale magnetic fluxes of fully-convective stars derived from spectropolarimetric data and Zeeman-Doppler Imaging (ZDI), as a 
function of mass (left panel) and mass xP rot 2 (right panel). Symbols are similar to those used in the mass-period diagram (see Fig.[TJ. For stars in the WM 
regime symbols corresponding to different epochs for a given star are connected by a vertical grey line. The yellow region represents the domain where 
bistability is observed and the orange one separates the two types of magnetic fields identified (see text). 




10" 1 

Mass (M Q ) 



Figure 3. Total magnetic fluxes of fully-convective stars in the saturated 
regime mea sured from unpolaris ed sp ectra of FeH lin e s. The values are 
taken from Rei ners et alj <2009t) a ndlReiners & Basril feOld) . whenever 
2MASS near infrared luminosities JCutri et alJl2003F and Hipparcos par- 
allaxes j ESA| 1 1997h are available to compute the stellar mass from the 
Delfos se et al mass-luminosity relation. Whenever spectropolari- 

metric data are available the properties of the magnetic topology are repre- 
sented as symbols described in Fig. [T] The magnetic field (y-axis) scale is 
the same as Fig. [2] The yellow region represents the mass domain where 
bistability is observed in spectropolarimetric data (see Fig. [2). 



tial term in the momentum equation. This weak field force balance 
requires small lengthscales. On the strong field branch, however, 
the Lorentz force relaxes the rotational constraint. A similar bifur- 
cation diagram, but based on the fact that magnetic buoyancy would 
be negligible close t o the d ynamo onset has been proposed for stars 
bv lWeiss & Tobias! feOOCh . 

The existence of a strong field dynamo regime has re- 
ceived support from theoretical and numerical studies (e.g., 
Childress & Sowardl fl973 : iFautrelle & Childress! I"l 9821 ; 1st. Pierrel 
19931) . More recently numerical simulations have supported the ex- 
istence of both branches in spherical geometry, i.e. both weak and 
strong field solutions were obtained depending on the initial condi- 
tions (Dormy & V. Morin, in prep. ). 
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Figure 4. Ant icipa ted bifurcation diagram for the geodynamo (adapted 
from [Roberts 1988). The magnetic field amplitude is plotted against 
the Rayleigh number. The bifurcation sequence is characterised by two 
branches, referred to as weak and strong field branches. The yellow and or- 
ange regions have the same meaning as in Fig. [2] Rs^ is the critical Rayleigh 
number for the onset of non-magnetic convection. The weak field regime 
sets in at Ra m , and the turning point associated with the runaway growth 
corresponds to Ra = Ra r . 



4 DISCUSSION 

We now speculate that the group of stars showing multipolar and 
time-varying magnetic topologies (WM) correspond to the weak 
field regime, whereas those with a steady dipole (SD) belong to the 
strong field branch. 

The usual control parameter in the weak vs strong field dy- 
namo scenario described above is the Rayleigh number, which mea- 
sures the energy input relative to forces opposing the motion. Mass 
can be used as a good proxy for the available energy flux in M 
dwarfs, Fig. [2^ can therefore be interpreted as a bifurcation dia- 
gram for the amplitude of the large scale magnetic field versus a 
control parameter measuring the energy input. In order to com- 
pare the driving of convection with the impeding effect of rota- 
tion, we can use MP^,, as a rough proxy for the Rayleigh num- 
ber (see Fig. [2})) based on rotation rather than diffusivities (e.g., 
IChristensen &Aubertll2006l) . 

Such an identification implies that (i) the strong field/weak 
field dichotomy only affects the large-scale component of the mag- 
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netic field, (ii) the field strength is compatible with a Lorentz- 
inertia force balance for stars featuring a WM magnetism, whereas 
a Lorentz-Coriolis balance prevails for stars in the SD group. It is 
difficult to quantify the range of control parameter over which both 
branches co-exist, we will therefore focus our discussion on the 
prevailing force balances in both regimes and their implications on 
the magnetic field. 



4.1 Large-scale dynamo bistability 

Different types of magnetism have previously been found to affect 
only the large-scale component of the magnetic field (measured 
with spectropolarimetry) and not the total magnetic flux (mea- 
sured with unpolarised spectroscopy). Indeed, the aforementioned 
spectropolarimetric survey has revealed that the large-scale mag- 
netic field of M dwarfs rapidly changes with stellar mass (both 
in geometry and field strength) close t o the fully-convective limit 
(Morin et al. 2008bl; lDonati et alj|2008l) , w hereas no change is vis - 
ible in total magnetic flux measurements jReiners & Basril l2007l) . 
As the large-scale component only represents a small fraction of 
the total flux, a change affecting the large-scale field alone can in- 
deed remain unnoticed in the measured values of the total field. 

The Rossby number in stars is much higher than in the Earth's 
interior and associated with a stronger driving. Therefore, whereas 
the geodynamo must act on comparatively large scales (because of 
the fairly moderate value of the magnetic Reynolds number) mo- 
tions in stellar interiors most likely generate fields on a variety 
of scales, this includes the possible c oexistence of a large scale 
dynamo with a small sc ale dynamo jVogler & Schusslen 120071 : 
ICattaneo" & Hughes 2009). Such co-existence could easily account 
for the difference in measurements provided by spectropolarimetry 
(Fig.0 and unpolarised spectra (Fig.[3j- 



4.2 Force balance and magnetic field strength 

In the strong field regime the Lorentz and Coriolis forces are of 
the same order of magnitude. The ratio of the two forces can be 
estimated by the Elsasser number: 



A : 



B 2 
ppj\Q. 



(1) 



where B is the magnetic field strength, p the mass density, [i the 
magnetic permeability, 77 the magnetic diffusivity and fl the rotation 
rate. 

This magnetostrophic force balance is valid for large spatial 
scales which are strongly affected by the Coriolis force, and does 
not apply to small spatial scales for which the inertial term is pre- 
dominant in the momentum equation. It is important to note here 
that the Elsasser number only provides a crude measurement for 
this force balance. To establish this measure, an equilibrium be- 
tween induction and diffusion is assumed in the induction equa- 
tion. In doing so, the typical length scale of the field and of the flow 
have to be considered equal. While this is a sensible approximation 
for a planetary dynamo working at moderate magnetic Reynolds 
number, it is a crude approximation for stellar interiors. More im- 
portantly, this force balance can only provide an order of magnitude 
estimate for the field strength. The magnetic energy on the strong 
field branch will obvious ly vary with th e amount of thermal energy 
available (see Fig.|4]and Roberts 1988). Let us nevertheless try to 
provide an estimate of the surface field corresponding to an Elsasser 
number of unity for M dwarfs. We simply take p = M+ /(|w/J t 3 ), 



and similarly to C09 we assume that the ratio between the mag- 
netic field inside the dynamo region and the surface value is equal 
to 3.5. An estimate for the turbulent magnetic diffusivity is also re- 
quired, crude values, which can be derived from sunspot or active 
regions decay time or from the formula 77 ~ u rms £ (where u rms is the 
turbulent velocity and £ is a typical length scale ) are in the range 
10 1 1 - 3 x 10 12 cm 2 s" 1 (e.g., iRudiger et alj201ll) . Let us introduce 
rj m f = 10" cm 2 s . With r\ oc u rms f, assu ming that u rms scales with 
L* ' according to mixing-length theory dVitenseil953l) , and i with 
the depth of the convective zone, we derive an estimate of the field 
strength at the stellar surface in the strong field regime: 



(2) 



Taking stellar radius an d luminosity for the s tellar m ass in the range 
0.08 - 0.35 M from IChabrier & Baraffd d 19971) main sequence 
models, we note that B s f is almost independent of mass in this 
range, and thus the main dependence is on the rotation period and 
the chosen reference magnetic diffusivity rj e . 

However, we do not find evidence for a dependence of the 
large-scale magnetic flux on rotation period among stars belonging 
to the strong field branch in the spectropolarimetric data. Depend- 
ing on the precise extent of the bistable domain, a factor of 2-3 in 
magnetic fluxes would be expected between the fastest and slow- 
est rotators of our sample. Such a moderate dependence might re- 
main undetected due to the dispersion (object-to-object variations) 
of measurements. Using the aforementioned estimate of rj a and the 
rotation periods of the stars in the spectropolarimetric sample, we 
find surface values in the strong field regime ranging from 2 to 
50 kG. Such estimates are compatible with the order of magnitude 
of the measured large-scale magnetic fluxes. It should be stressed 
however that this is not conclusive as the weak field branch is only 
a factor of ten smaller in magnitude. 

The weak field vs strong field scenario can however receive 
additional support by considering the ratio of field strengths on the 
weak and strong field regimes (i.e. the amplitude of the gap be- 
tween both branches). In the case of the Earth's core, inertia and 
viscous terms become significant at a similar length scale (because 
Ro 2 E~' ~ 0(1)). In stellar interiors, however, the Rossby num- 
ber is usually much larger than in the Earth's core (by a factor 
of at least 10 4 ). A weak field branch would therefore naturally re- 
sult from a balance between Lorentz and inertial forces through the 
Reynolds stresses (i.e. the field strength on the weak field branch 
should approach equipartition between kinetic and magnetic ener- 
gies). Therefore, the ratio of field strengths on the weak and strong 
field regimes - corresponding to a different force balance - is ex- 
pected to depend on the ratio between inertia and Coriolis forces. 
We can thus estimate: 



Ro 



1/2 



(3) 



where Ro is the Rossby number. MlOa derived empirical Rossby 
numbers of the order of 10~ 2 for stars in the bistable domain, im- 
plying a ratio B s f/B w f of the order of 10, which is indeed the typical 
ratio of large-scale magnetic fields measured between the WM and 
SD groups of stars (see Fig.|2]and sec. [2J. 

We shall now discuss an apparent caveat of the proposed sce- 
nario. As noted in section [2] all the stars considered in our sample 
belong to the so-c alled saturated regime dKiraga & StepierJl2007l: 
iReiners et alj2009l) . This means that the strong variation of the field 
strength with P Iot observed for lower rotation rate has "saturated" 
and the amplitude of the field seems to be independent of the rota- 
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tion rate of the star. Hence, there is an apparent contradiction with 
the possibility of a strong field branch, on which the magnetic field 
depends on rotation rate as B s f oc fl 1/2 . The first important point 
is that the B sf oc SI 1 ' 2 (derived from A ~ 1) should apply here 
to the large scale field alone, which is only a fraction of the total 
magnetic field of the stars (between 15 and 30 %). If a small scale 
dynamo operates, it does not need to follow the same dependency. 
Besides, the slope of this flat portion of the rotation-magnetic field 
relation (either of the overall magnetic flux based on unpolarised 
spectroscopy or of its proxy the relative X-ray luminosity L x /L bol ) 
is poorly constrained in the fully-convective regime and is in fact 
compatible with an 1 ' 2 dependence. Evidence for such a depen- 
dence of the large-scale magnetic field on fl 1 ' 2 would strongly sup- 
port the proposed scenario. 



5 CONCLUSION 

In the present letter, we compare the bistability predicted for the 
geodynamo with the latest results on spectropolarimetric observa- 
tions of fully convective main-sequence stars (MlOa). We show that 
the weak vs strong field dynamo bistability is a promising frame- 
work to explain the coexistence of two different types of large- 
scale magnetism in very low mass stars. The order of magnitude 
of the observed magnetic field in stars hosting a strong dipolar 
field (SD), and more conclusively the typical ratio of large-scale 
magnetic fields measured in the WM and SD groups of stars are 
compatible with theoretical expectations. We argue that the weak 
dependency of the magnetic field on stellar rotation predicted for 
stars in the strong-field regime cannot be ruled out by existing data 
and should be further investigated. We do not make any prediction 
on the extent of the bistable domain in terms of stellar parameters 
mass and rotation period, this issue shall be investigated by further 
theoretical work, and by surveys of activity and magnetism in the 
ultracool dwarf regime. 

A dynamo bistability offers the possibility of hysteretic be- 
haviour. Hence the magnetic properties of a given object depend 
not only on its present stellar parameters but also on their past evo- 
lution. For instance, for young objects episodes of strong accretion 
can significantly modify their structure a nd hence the conv ective 
energy available to sustain dynamo action dBaraffe et alj2009T) ; ini- 
tial differences in rotation periods of young stars could also play a 
role. Because stellar magnetic fields are central in most physical 
processes that control the evolution of mass and rotation of young 
stars (in par ticular accretion-ejections processes and star-disc cou- 
pling, e.g., lBouviej r2009: Gre gory et al] [20 10). the confirmation 
of stellar dynamo bistability could have a huge impact on our un- 
derstanding of formation and evolution of low mass stars. 



ACKNOWLEDGEMENTS 

It is a pleasure to thank Katia Ferriere for helpful discussions and 
comments on this manuscript. We acknowledge the constructive 
comments of an anonymous referee. 



REFERENCES 

Baraffe I., Chabrier G., Gallardo J., 2009, ApJ, 702, L27 
Bouvier J., 2009, in C. Neiner & J.-P. Zahn ed., EAS Publications 
Series Vol. 39 of EAS Publications Series, The role of magnetic 



fields in governing the angular momentum evolution of solar- 
type stars, pp 199-209 
Browning M. K., 2008, ApJ, 676, 1262 
Cattaneo F, Hughes D. W., 2009, MNRAS, 395, L48 
Chabrier G., Baraffe I., 1997, A&A, 327, 1039 
Chabrier G., Kuker M., 2006, A&A, 446, 1027 
Chandrasekhar S., 1961, Hydrodynamic and hydromagnetic sta- 
bility. International Series of Monographs on Physics, Oxford: 
Clarendon, 1961 
Charbonneau P., MacGregor K. B., 1997, ApJ, 486, 502 
Childress S., Soward A. M., 1972, Phys. Rev. Letters, 29, 837 
Christensen U. R., Aubert J., 2006, GeoJI, 166, 97 
Christensen U. R., Holzwarth V., Reiners A., 2009, Nature, 457, 
167 

Cutri R. M. et al., 2003, 2MASS All Sky Cata- 
log of point sources.. The IRS A 2MASS All-Sky 
Point Source Catalog, NASA/IPAC Infrared Science 
Archive. http://irsa.ipac.caltech.edu/applications/Gator/ 

Delfosse X. et al., 2000, A&A, 364, 217 

Donati J.-F. et al., 2003, MNRAS, 345, 1 145 

Donati J.-F, Brown S. F, 1997, A&A, 326, 1135 

Donati J.-F, Forveille T, Cameron A. C, Barnes J. R., Delfosse 
X., Jardine M. M., Valenti J. A., 2006, Science, 311, 633 

Donati J.-F. et al., 2006, MNRAS, 370, 629 

Donati J.-F. et al., 2008, MNRAS, 390, 545 

Eltayeb I. A., Roberts P. H, 1970, ApJ, 162, 699 

ESA 1997, VizieR Online Data Catalog, 1239, 

Fautrelle Y, Childress S., 1982, Geophys. Astrophys. Fluid Dyn., 
22, 235 

Featherstone N. A., Browning M. K., Brun A. S., Toomre J., 2009, 
ApJ, 705, 1000 

Goudard L., Dormy E., 2008, Europhysics Letters, 83, 59001 
Gregory S. G., Jardine M., Gray C. G., Donati J., 2010, Reports 

on Progress in Physics, 73, 126901 
Kiraga M., Stepien K., 2007, Acta Astronomica, 57, 149 
Morin J., et al., 2010, arXiv, arXiv: 1009.2589 
Morin J., Donati J., Petit P., Delfosse X., Forveille T, Jardine 

M. M., 2010, MNRAS, 407, 2269 
Morin J. et al., 2008, MNRAS, 384, 77 
Morin J. et al., 2008, MNRAS, 390, 567 

Phan-Bao N., Lim J., Donati J.-F, Johns-Krull C. M., Martin 

E. L., 2009, ApJ, 704, 1721 
Reiners A., Basri G., 2006, ApJ, 644, 497 
Reiners A., Basri G., 2007, ApJ, 656, 1 121 
Reiners A., Basri G., 2010, ApJ, 710, 924 
Reiners A., 2010, arXiv, arXiv: 1012. 1183 
Reiners A., Basri G., Browning M., 2009, ApJ, 692, 538 
Roberts P. H., 1978, Rotating fluids in geophysics 
Roberts P. H., 1988, Geophys. Astrophys. Fluid Dyn., 44, 3 
Rudiger G., Kitchatinov L. L., Brandenburg A., 2011, Sol. Phys., 

269,3 

Saar S. H., 1988, ApJ, 324, 441 
Semel M., 1989, A&A, 225, 456 

Simitev R. D., Busse F. H., 2009, Europhysics Letters, 85, 19001 
Soward A. M., 1979, Phys. Earth Planet. Int., 20, 134 
St. Pierre M. G., 1993, in M. R. E. Proctor, P. C. Matthews, & 
A. M. Rucklidge ed., Solar and Planetary Dynamos The strong 
field branch of the Childress-Soward dynamo., pp 295-302 
Stevenson D. J., 1979, Geophys. Astrophys. Fluid Dyn., 12, 139 
Vitense E., 1953, ZAp, 32, 135 
Vogler A., Schussler M., 2007, A&A, 465, L43 
Weiss N. O., Tobias S. M., 2000, Space Sci. Rev., 94, 99 



